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Ph ' Using Husinii function approach, we study the "quantum phase space" 

of a harmonic oscillator interacting with a plane monochromatic wave. We 
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show that in the regime of weak chaos, the quantum system has the same 



^ ' symmetry as the classical system. Analytical results agree with the results of 

^" ' numerical calculations. 
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- ■ It is known that the phase space of a classical harmonic oscillator weakly interacting with 

'^i a plane monochromatic wave possesses an interesting symmetry. (See, for example,a and 

-(— > 

S . references therein.) In the case of exact resonance, /icu = fi (/i = 1, 2, . . .), between the wave 

^. 

O^. (with the frequency Vt) and the harmonic oscillator (with the oscillation frequency cj), and 



under the condition e ^ 1 (where e is a dimensionless perturbation parameter), the classical 
phase space consists of an infinite number of resonant cells with the symmetry 2/i. An 
example of a corresponding phase space with /x = 4 is shown in Fig. 1. At the center of each 
cell there is an elliptic stable point. The particles move in the phase space around this point 
along the closed trajectories. The cells are separated from each other by the separatrices 
which are schematically shown in Fig. 1 by dashed lines. These separatrices form in the 
phase space an unlimited net. The net is covered by the stochastic layers forming the infinite 
stochastic web. When the perturbation parameter, e, is small the web width is exponentially 
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thinEl. However, if the particle is initially placed inside a stochastic region, it can travel 
throughout the web and gain energy, even for an arbitrarily small perturbation parameter, 
e. The existence of the crystalline and quasi- crystalline symmetries of the classical phase 
space, and stochastic web differ significantly this system from classical nonlinear systems 
with chaotic behavior .lil These interesting properties of the classical harmonic oscillator in a 
monochromatic wave motivated our studies of the corresponding properties in the quantum 
system. 

The quantum harmonic oscillator interacting with a monochromatic wave is described 
by the Hamiltonian, 

H = ^+ +-cos{kx-nt) = Ho + V{x,t), (1) 

where ^o is the Hamiltonian of the harmonic oscillator; V(a;, t), is the interaction Hamil- 
tonian; ejk and k are, respectively, the amplitude and the wave vector of the wave, x and 
f) are the coordinate and the momentum operators of the particle, m is the mass of the 
particle. The Hamiltonian (1) appears, for example, when analyzing the stability of an ion 
in a linear ion trap in the field of two laser beams with close frequencies^. The dynamics 
of the quantum system described by the Hamiltonian (|T]) is controlled by four parameters: 
the resonance number, /x; the detuning from the exact resonance: ^lo = fiuJ — Q; the di- 
mensionless perturbation parameter: e = ek/muj'^] and the dimensionless Planck constant: 
^0 = {hk"^) / (muj) . (For the system considered inB, the Lamb-Dicke parameter, 77, is related 
to fiQ-. TiQ = ^tf). Influence of these parameters on the dynamics of quantum system was in 
detail considered in Refs.B^ 

In this paper, we study the "quantum phase space" of a system described by the Hamil- 
tonian (|l]) for the case of exact resonance, 5uj = {] (when the number of the resonance cells 
is infinite) and under the condition: e <C 1 (when the chaotic layers, covering the separatrix 
net, are exponentially thin). We investigate the structure of the "quantum phase space". 
We show that quantum system possesses the same symmetry as the classical system. The 
structure of the quantum system with the Hamiltonian (0) is characterized by the Floquet 
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FIG. 1. The classical phase space for a harmonic oscillator in a monochromatic wave, for: 5uj = 
(exact resonance), ^ = 4, e = 0.05. 



states (or quasienergy states) found in Refs.ua. In order to build the phase space for the 
quantum system we use the Husimi functions of the Floquet states. The Husimi function for 
the wave function ipi^x, t) is defined as the projection of ip{x, t) on the coherent wave packet, 
xix; X, P), with the maximum at the point (X, P)Q, 



$(X, P, t) = ^\< tp{x, t)\x{x; X, P) > |2. 



(2) 



The Husimi function, $(X, P, t), defines the probability of finding a quantum particle char- 
acterized by the wave function tp{x,t) at the point {X ,P) of the "quantum phase space". 



The cross-sections of the Husimi function are the hnes of equal probabihty of finding the 
quantum particle. Below, these lines for Husimi functions are compared with the trajectories 
in classical phase space. 

Namely, we analyze the structure of the Husimi functions of the quasienergy states and 
compare them with the structure of the classical phase space. First, we present some general 
formulas which will be used to investigate the system described by the Hamiltonian (|l]). It is 
convenient to decompose the coherent state, xi^] ^y P)y i^ito the complete set of harmonic 
oscillator eigenstates. 



X{x; X, P) = exp ( ) }^ ^ _ > m(a^), 



4^0 / m=o ^{2ho)^m\ 



(3) 



where ipn{x) is the n-th eigenfunction of the harmonic oscillator with the Hamiltonian ifo- 
In Eq. (^) we used the dimensionless coordinate (X = kx) and the dimensionless momentum 
(P = pk/muj). We use the same basis to represent of the wave function, ip{x,t), 

1 



Ipix, t) = Y^ Cn{t)i}n{x) exp 



n=0 



-iiot n + 



(4) 



The structure of the Husimi function (^) is completely defined by the coefficients Cn{t), 



$(X, P, t) 
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5:Q(t)i^=I=4=exp(^mu;t) 
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It is convenient to use cylindric coordinates. 



X = rcos(/9, P = rsin(/?, 



(6) 



where r = yX^ + P^ and 99 = arctg(P/X). In these variables, the Husimi function (^ is. 



$(r, (^, t) 



exp 
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J2 Q(^) n ^ =exp(im^t) 

m=o J{2ho)"'m\ 



(7) 



Since the perturbation, V{x,t), in (1) is periodic in time, one can use Floquet theory 
and write the solution of the non-stationary Schrodinger equation as. 



ipg{x,t) = exp{-iEgt/h)Ug{x,t), 



(8) 



where Uq{x,t) = Uq{x,t + T) is a time-periodic function whose period is T = 27r/r2. The 
index q labels the quasienergy (QE) states. It is convenient to use the complete set of 
harmonic oscillator eigenf unctions to represent the function Uq{x,t), 

oo 

Uq{x,t) = Y.cmMx), (9) 

n=0 

where the expansion coefficients, C^{t) = C^(t + T), are time-periodic functions. Using Eqs. 
(0)"(S) ^^ '^^^ rewrite the Husimi function of the QE state as, 

, ,.2 . 2 
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y2 C^ I — exp ( 2mms 

m=o J{2ho)"'m\ 



(10) 



where s = 0, 1, 2,.... The QE states of the monochromatically perturbed harmonic oscilla- 
tor were studied in detail in a series of papergai^, using degenerate resonance perturbation 
theory for the Floquet states. In particular, the quantum regimes corresponding to regular 
motion and to the case of weak chaos in the classical phase space were investigated. As it was 
shown in Ref.a, the Hilbert space of the quantum system breaks up to some approximation 
into the dynamically independent regions — quantum resonance cells, each of them with 
its own set of QE states. In the zeroth order (resonance) approximation, the QE functions 
and the QE spectrum of each cell are almost independent. Near the top and bottom of the 
QE spectrum of an individual cell, the QE states are the states of an effective harmonic 
oscillator. The QE levels are equally spaced, with a separation huj between the levels. The 
frequency, uj, in the quasiclassical limit coincides with the frequency of small oscillations 
near the center of the resonance in phase space. In this paper, we consider only two extreme 
QE states of an individual cell - extreme upper and extreme lower QE states, called the 
"QE ground states". Thus, each quantum resonance cell has two QE ground states. The 
QE functions, C^, and the QE levels. Eg, of the ground states of each individual cell are 
connected by the relationsQ, 

Eq ~^ ~Eg, C^^ -^ (— 1)"'C|J^, (11) 

associated with the transformation: x — »• —x in Eqs. (§),(§). 



Below we analyze the structure of the Husimi functions of the QE ground states. The 
upper QE function is a Gaussian wave packet, 

where F is the normalization factor, and Ue is the position of the maximum of the QE wave 
packet in the Hilbert space which corresponds to the quantized radius of the elliptic stable 
point: Te = \/2neho (see Ref.a). The width of the wave packet, Oe, in Eq. ( |12|) was defined 
in Ref-Q in the form, 

were the function g^{n) is expressed in terms of the matrix element: gfj,{n) =< 
ipnl cos{kx)\tpn+^i >■ In the quasiclassical region of parameters, Eq. (|^) can be expressed 
through the half width in action of the classical resonance cell. A/ (expression for the value 
of A/ see for example in Ref .0) , 
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1/4 / A r \ 1/2 



(14) 
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where the prime indicates differentiation with respect to the argument. For example, for 
/i = 1 we have: Oe = ''"e/{^o[(^e)^ — l]}^''^- The boundaries of the quantum cells are given by 
the zeroes of the function g^{n). As was shown in Ref.H, the function g^{n) is proportional 



to the Bessel function, J^, of order /i: ^'^(tt,) ~ J^{\^2nho). So, the number of levels in the 
individual cell is proportional to ho. Thus, the ratio: (the packet's width in n)/(the cell's 
width in n) is proportional to v^o? s^nd in the quasiclassical limit the relative width of the 
QE ground state tends to zero. 

The Husimi representation allows one to construct the QE eigenstates in the quantum 
phase space. The simplest case is the Husimi function of a single harmonic oscillator state: 
Cn = (^n,no) which due to Eq. (j^) has the form, 

, , exp(-77^) r^"o 
<l>("")(r, ^, t) = ^ ^'''\ ,\ -. (15) 



This expression has its maximum at tq = ■\/2no^o- The definite value of uq corresponds to 
the definite value of the action: Jq = hoUQ. Due to the fundamental uncertainty relation, 
the phase, ip, of this state is indefinite. The Husimi function is independent of the phase, ip, 
and looks like a round hump. 

In agreement with Eqs. (10) and (12), the Husimi function of the ground QE state is, 

2 
(m — JT-e)^' 
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(16) 



Only Am ~ a^ terms with \m — ne\ < 2ae effectively contribute to the sum on the right-hand 
side of Eq. (16), and one can neglect all other terms. Then, Eq. (|l^) becomes. 



where we assumed: ng 3> 1, so that. 



Am ^ngm(v+2^) 

E ;,_. , exp 



n=-Am J{2hone)'^ 
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2al 



(17) 



(rie + 'Ti)! ~ nglrig 



(18) 



The double sum in Eq. (O) can be rewritten as. 



Am ^„+mgi(n-m)((/3+-^J 



n,m=-Am y(2/lo^e)"^ 



exp 






2Am 
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E I /H^— r l ^^P(-ib) E e^'^^^+'Dexp 
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j=-2Am \y2Tione) \ T:".e/ fc=-2Am \ """e , 

where j = n + m, k = n — m. Thus, by using the approximation (|l^) we find that the 
Husimi function of the extreme QE state, can be factored. 



$e(r, ^,s) =7(r)^((/?,s). 



(19) 



InEq. (19), 
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We now find the coordinates of maxima of $e(^, f)- Suppose that each maximum of the 
Husimi function corresponds to the stable elliptic point at the center of a resonance cell. 
Maximum of ^q{r, if) in r is defined from the equation, 

d e 2^0 r^'^Hrl' 



d 
dr 



dr 27r 2/^0"' ^e' 



,11©'°-^"^)^ 
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When r = Te, both sums in Eq. (^) are zero: in the first term, the derivative is equal to 
zero as follows from Eq. ([TsD ; in the second term, the sum is equal to zero, and the value r^ 
can be considered as the radius of the center of the quantum resonance cell in the quantum 
phase space. 

We now find the maxima of ^(^9, s). It is convenient to present this function in the form, 

e(^) = 1 + 2 Y.^ cos(/.m^) exp (^-^^ j , (23) 

where we took into account that in the resonance approximation the particle can populate 
only states with the numbers: k = fxin (see Ref.a). All terms in the sum on the right-hand 
side of Eq. ( p3D decrease in absolute values as m increases. Then, the extrema of the 
function ^{^p) is defined by the the extrema of the term with m = 1. When /x = 1 there is 
one maximum at y? = 0; when fi = 2 there are two maxima ai ip = and ip = tt. In general 
case the function C,{f) has /i maxima. 

The extreme lower QE function is related to the extreme upper one by the transformation 
(p!T|), which is convenient to rewrite in the form: C^^ -^ exp(— ?7rm)C^„. The function 
C,iower{^) of the lower ground QE state is, 

iloweri,^) = 1 + 2 Yl COS [(/!(/) - T[)m] CXp I — ^- 1 . (24) 



The maxima of the function ^(v?) in Eq. (|23D correspond to minima of ^loweri.'P) in Eq- (^4|), 
and vice versa. Thus, for /x = 1 the function ^iower{'^) has a maximum at (y9 = vr; at yU = 2 
there are two maxima at </? = ±7r/2 and so on. In general, the Husimi functions of the two 
QE ground states have 2/i maxima with the radius r = r^. Each maximum is situated at 
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FIG. 2. Contour plots of the Husimi functions for the exact resonance case, with the resonance 
number ^ = A; ho = 0.12; e = 0.002. 



the center of a quantum resonance cell, so that the quantum phase space has the same 
symmetry as the classical phase space. For /x = 4, the symmetry of the Husimi function, 
shown in Fig. 2, is the same as the symmetry of the classical phase space in Fig. 1. A 
similar result was demonstrated numerically in Ref.Q for fi = 1. As one can see from Fig. 
2, in agreement with Eq. (11), the quantum phase space is symmetric with respect to 
the substitution: X -^ —X. However, there is no exact symmetry with respect to the 
transformation: P ^ —P. The reason is presumably related to our approximation (p!8D 



which leads to separation of variables in Eq. (|T^). This approximation is more valid for the 
"quasiclassical cells" with n ^ 1 than for "quantum cells" , for which the value of n is not 
large. One can see from Fig. 2, that "quasiclassical cells" are more symmetrical than the 
"quantum cells" , and the structure of the "quasiclassical cells" is close to the structure of 
the classical cells shown in Fig. 1. This symmetry of the quantum phase space differs this 
system from quantum chaotic systems with critical threshold to global chaos.l3 

In summary, the correspondence between the symmetry of the Husimi functions of the 
QE ground states and the symmetry of the classical phase space has been demonstrated for 
a degenerate system both analytically and numerically. 
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